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Lattice Study of the Massive Schwinger Model with a 6 term under 
Liischer's "Admissibility" condition* 

H. Fukaya and T. Onogi'*^ 

'^Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

We present a numerical study of the massive two-flavor QED in two dimensions with the gauge action proposed 
by Liischer, which allows only "admissible" gauge fields. We find that the admissibility condition does not allow 
any topology changes by the local updation in Hybrid Monte Carlo algorithm so that the configurations in each 
topological sector can be generated separately. By developing a new method to sum over different topological 
sectors, we investigate 6 vacuum effects. Combining with domain-wall fermion action, we obtain the fermion 
mass dependence and 9 dependence of the meson masses, which are consistent with the analytic results by mass 
perturbation in the continuum theory. 



1. Introduction 

The studies of chiral symmetries on the lattice 
have made much progress recently. There are two 
keys in these developments; (1) Ginsparg- Wilson 
relation [1]: 



751? + Dj5 = aDj^D, 



(1) 



and (2) Liischer's "admissibility" condition [2,3]: 



<e for all x,fi,,^. 



(2) 



where a denotes the lattice spacing, U^^{x) de- 
notes a plaquette and e is a fixed constant. 

The first key, the Ginsparg- Wilson relation, 
gives a redefinition of chiral symmetries without 
fermion doublers at the classical level. The sec- 
ond key, the "admissibility" condition Eq. (2), 
makes gauge fields smooth. Under this condi- 
tion the gauge fields are classified by a topological 
charge which corresponds to that of the contin- 
uum theory, giving a well-defined chiral anomaly 
at the quantum level. 

In two-dimensional QED, the topological 
charge is defined as 

Q = 2^E.,^,.[-f£p.lnC/;.W] (integer) 
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While numerical studies on Ginsparg- Wilson 
fermion have been performed, no study of 
Liischer's admissible gauge fields has been done 
before. 

We present a numerical study of two- 
dimensional QED with Liischer's gauge action 
which satisfies admissibility condition automat- 
ically [4]. We show that topological structure is 
realized on the lattice, 9 dependence of observ- 
ables can be evaluated by our new method. As 
a result we obtain the mass of isotriplet meson 
which is consistent with the continuum theory. 
This work may offer a new approach to investi- 
gate "topology" in the lattice gauge theories. 



2. Lattice Simulations 

We take Liischer's gauge action combined with 
the domain-wall fermion action; 



PSc + ^ Sf, 
fiavor=l 



So - 



St 



-{l-RcUPix))/e 



if admissible 
oo otherwise 

the domain- wall fermion action 

including Pauli-Villars fields. 



(4) 



(5) 



(6) 
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where (3 = The Liischer's condition is au- 

tomatically satisfied with this action. 

We use the Hybrid Monte Carlo (HMC) 
method with pseudo-fermions and conjugate gra- 
dient (CG) algorithm for the inversion of the 
fermion matrix. 50 molecular dynamics steps 
with step-size At = 0.02 are performed for each 
trajectory and configurations are updated per 10 
trajectories. We also check the admissibility at 
the Metropolis test. The simulations are carried 
out on a 16 X 16 X 6 lattice with /3 = 0.5 and 
e = 1.0. Fermion mass m is taken to be 0.1,0.2,0.3 
and 0.4. 

3. Topological Sectors and Reweighting 

We set the initial link variables as follows. 



2TrNi , , 
exp —Ox,Ly \ , 



Ur\x,y) = expj^.}, (7) 

where L = 16 is the lattice size. This configura- 
tion gives constant background electric fields over 
the torus and minimizes the gauge action in the 
sector with topological charge A''. As Fig.l shows, 
Liischer's action generates configurations without 
any topology changes. 

Thus admissible gauge fields are separated into 
topological sectors and configurations are gener- 
ated in each sector with Liischer's action. The 
next task is to sum up the expectation values in 
each sector with correct weights. 



(0)f" 



E+oo 



//3,m 



(8) 



where 



i?,"(/3,m) 



Zo{l3,m) 

2xpcctatioi 

ZN{l3,m) : generating functional in each sector. 



(O)^ m • expectation value in each sector 



We developed a new method to calculate the 
reweighting factor m) [4]. By separating 

R^{(3,m) into three factors as 

ii^(/3,m)=exp(-/35^„,J 
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Figure 1. The comparison of the evolution of the 
topological charge with Wilson's gauge action and 
Liischer's gauge action for the same lattice spac- 
ings determined from the string tension. 
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(9) 



where SQ^-^j^, in the first factor denotes the mini- 
mum of the gauge action in N sector, the second 
factor is the contribution from free fermion deter- 
minants given by the integral of Polyakov loops , 
and the third one is a /? integral of the expectation 
value of the gauge action. 

Once the reweighting factor is known, we can 
evaluate theta dependence of observables as in 
Eq. (8). 
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4. Meson Masses 

We calculated the isotriplet meson propagators 

at various 9 as 



^(7r(a;,y)7r(0,0))/«H = 



y 

4 

^ e^^«^(7r(x,y)7r(0,0))^„i?^(/3,m). (10) 

JV=-4 y 

Here we have ignored |A^| > 4 sectors since they 
only give contributions less than 1.2 % of zero sec- 
tor. Fitting them into hyperbolic cosine function, 
we evaluated the meson mass. 

As Figs. 2 and 3 show, the fermion mass and 9 
dependence of the pion mass are consistent with 
that of continuum theory at small 6 [5-7] ; 



rriTT oc (to cos 2) ^ ■ 



(11) 



Moreover, it is remarkable the residual mass in 
the chiral limit is consistent with zero within sta- 
tistical error even in this strong coupling regime; 



m^(m ^ 0) = -0.057 ± 0.060. 



(12) 
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Figure 2. The fermion mass dependence of the 
pion mass at ^ = 0. 
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5. Summary 

We performed a numerical simulation of 2 di- 
mensional theory with Liischer's gauge action 
combined with domain-wall fermion action. We 
also proposed a new method for summing over 
different topological sectors. 

Applying Liischer's condition to numerical sim- 
ulations we find that the lattice gauge fields are 
separated into topological sectors. Our reweight- 
ing method enables the summation over topolog- 
ical sectors. We obtained the results which are 
qualitatively consistent with that in the contin- 
uum theory. It is also notable that the resid- 
ual mass of the domain-wall fermion in the chiral 
limit is consistent with zero even in the strong 
coupling regime. 

As future works, we plan to perform more pre- 
cise and quantitative studies of 2 dimensional 
QED and comparison to the results in conven- 
tional lattice formalism [8]. Studying how to ex- 
tend our approach to 4 dimensional QCD is also 
important. 
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Figure 3. 9 dependence of the pion mass at m = 
0.2. 
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